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Abstract 

Self-consistent account of the most simple non-gauge vector fields leads to a broad spectrum 
of regular scenarios of temporal evolution of the Universe completely within the frames of the 
Einstein's General relativity. The longitudinal non-gauge vector field is "the missing link in the 
chain", displaying the repulsive elasticity and allowing the macroscopic description of the main 
features of the Universe evolution. The singular Big Bang turns into a regular inflation-like state 
of maximum compression with the further accelerated expansion at late times. The parametric 
freedom of the theory allows to forget the troubles of fine tuning. In the most interesting cases the 
analytical solutions of the Einstein's equations are found. 
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I. INTRODUCTION 



Regular scenarios of the Universe evolution, driven by the non-gauge vector fields together 
with the ordinary matter, are considered within the frames of Einstein's theory of general 
relativity. 

From the standpoint of general relativity the matter curves the space-time, giving rise 
to mutual attraction between the bodies. However, according to modern observations, the 
Universe is expanding as a whole, despite the gravitational attraction between material ob- 
jects. The expanding solution of the Einstein's equations due to the cosmological constant 
belongs to De Sitter The expanding solutions of the Einstein equations without the 
cosmological constant (Friedman- Robertson- Walker (FRW)Q) inevitably contained the sin- 
gularity. The singularity works as a cover for the unknown hidden origin of expansion of 
the Universe, containing only mutually attracting material objects. For a long time the 
singularity is considered as a general property of the Universe. The singular point, referred 
to as "Big Bang", is commonly accepted as the "date of birth" of the Universe. Discovery 
of the accelerated expansion of the Universe shows that the source of acceleration continues 
to exist for a long time after the Big Bang. Naturally, the fact of accelerated expansion 
gave rise to the assumption that the physical vacuum is not just the absence of the ordinary 
matter. The existence of the so called " dark energy" and " dark matter" , as the unknown 
source of the Universe expansion, is widely discussed in modern literature j^]. 

Among numerous attempts to guess the riddle of accelerated expansion, from my point 
of view, the most attractive one is the macroscopic description of the Universe expansion, 
driven by vector fields. Utilization of vector fields in general relativity shows undoubtable ad- 
vantages in comparison with scalar fields and with multiplets of scalar fields. The equations 
appear to be more simple, while their solutions are more general. The solutions have addi- 
tional parametric freedom, allowing to forget the fine-tuning problem [4j. However, starting 
from the pioneer paper by Dolgov[5[, people considered mostly gauge vector fields 6]- 10] in 
applications to the dark sector. 

Historically in flat space-time the divergence of the vector field was artificially set to zero 



ll| in applications to elementary particles. This restriction allowed to avoid the difficulty 



of negative contribution to the energy. In the electromagnetic theory it is referred to as 
Lorentz gauge. The gauge invariance takes place only for massless fields. The Lorentz gauge 
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does not allow to utilize all the advantages of vector fields. In general relativity (in curved 
space-time) the energy is not a scalar, and its sign is not invariant against the arbitrary 
coordinate transformations. Considering the vector fields in general relativity, it is worth 
rejecting the gauge restriction, using instead a more weak condition of regularity. Step by 
step, people are now getting rid of the Lorentz gauge restrictions [121]- [15]. 

Vector fields in general relativity form a three-parametric variety 4|. The analysis of 
the vector fields in the background of the De Sitter metric, including those with a nonzero 
covariant divergence, is performed in 16|. It came out that the energy-momentum tensor 
of the most simple zero-mass vector field enters the Einstein equations as an additive to the 
cosmological constant. Its back reaction affects the Hubble constant - the rate of expansion. 
The zero-mass vector field does not vanish in the process of expansion. On the contrary, 
massive fields vanish with time. Though their amplitude is falling down, the longitudinal 
massive fields make the expansion accelerated la ]. In other words, the longitudinal massive 
fields with a nonzero covariant divergence display the repulsive elasticity. It is worth analyz- 
ing the possible scenarios of the Universe evolution under the joint action of repulsing vector 
fields and attracting ordinary matter. Dynamical competition of repulsing and attracting 
forces results in the variety of regular scenarios of the Universe evolution. 

The paper is organized as follows. General properties of vector fields in the background 
of an arbitrary metric are presented in Section II. In Section III we present and analyze the 
Einstein equations, describing the evolution of the Universe under the action of the most 
simple non-gauge vector fields. The zero-mass vector field is responsible for either contrac- 
tion, or expansion, at a constant rate. The massive vector fields grow in the process of 
contraction and vanish with expansion. They are responsible for the regular transition from 
contraction to expansion, thus demonstrating the repulsive elasticity and ability to resist 
compression. In the limits of large and small masses the regular solutions are analyzed ana- 
lytically. This analysis facilitates considering the regular scenarios of the Universe evolution 
under the joint action of repulsing vector fields and attracting ordinary matter (Section IV). 
In most interesting cases the regular solutions are found analytically. The existence of four 
arbitrary dimensionless parameters allow to forget the fine tuning problem. In the limit of 
no vector fields the regular analytical solutions reproduce the Friedman-Robertson- Walker 
0] fine tuned singular cosmology. 

The results are discussed in Section V. 
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II. VECTOR FIELD IN GENERAL RELATIVITY 



In general relativity, the Lagrangian of a vector field <pi consists of the scalar bilinear 
combinations of its covariant derivatives and a scalar potential V{(f) K (f)K)- A bilinear combi- 
nation of the covariant derivatives is a 4-index tensor Sjklm — ( t ) i;K ( t ) L;M- The most general 
form of the scalar S, formed via contractions of Siklm, is S = (ag IK g LM + bg IL g KM + 
cg IM ' 9 )Siklmi where a, b, and c are arbitrary constants. The general form of the La- 
grangian of a vector field <pi is 

L = a(0^) 2 + 60f A/ 0f + c0f A/ 0f/ - V{<j> M <l> M ). (1) 

The classification of vector fields <pi is most convenient in terms of the symmetric Gjk = 
k(4>i-,K + ^K;i) an d antisymmetric F IK = \{4>i-,k — 4>K;i) parts of the covariant derivatives. 
The Lagrangian (CQ) gets the form 

L = a{G%f + {b + c)G L M G™ + (b - c)F L M F™ - V(^ M ). (2) 

The bilinear combination of antisymmetric derivatives Fj^F 1 ^ is the same as in electrody- 
namics. It becomes clear in the common notations Aj = 4>i/2, Fik = Aj-k — Ak-j. 

The terms with symmetric covariant derivatives deserve special attention. In applications 
;o elementary particles in the flat space-time the divergence is artificially set to zero 
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This restriction allows to avoid the difficulty of negative contribution to the energy. In the 
electromagnetic theory it is referred to as Lorentz gauge. However, in general relativity 
(in curved space-time) the energy is not a scalar, and its sign is not invariant against the 
arbitrary coordinate transformations. Considering the vector fields in general relativity, it 
is worth getting rid of the restriction ([3]), using instead a more weak condition of regularity. 
The covariant field equations 

^ K K;I + H'l K K + = -V'h (4) 

and the energy-momentum tensor 

T IK = -g IK L + 2V'0 7 K + 2ag IK ( ( f ) fl I <p L ) ;L + 2(6 + c)[(GW L ) ;L - G L K F lL - G\F KL \ 
+2(6 - c^F^Flk - F L K . L (pj - F\ L ct> K ) 

(5) 



describe the behavior of vector fields in the background of any arbitrary given metric gix 
3- Here V = f&^l. 

If the back reaction of the field on the curvature of space-time is essential, then the metric 



obeys the Einstein equations 



Rik - ^9ikR + J^9ik = xT IK (6) 



with added to Tjk- Here A and x are the cosmological and gravitational constants, 
respectively. With account of back reaction the field equations (TjJ are not independent. 
They follow from the Einstein equations ([6]) with Tjk © due to the Bianchi identities. The 
field equations (j3J) are linear with respect to <ft if the vector field is small, and the terms with 
the second and higher derivatives of the potential V (0m</> M ) can be omitted. 

III. VECTOR FIELDS IN COSMOLOGY 

Today it is generally accepted that among the staff of the Universe only 4.5% comes 



from the ordinary matter 18[. It is reasonable to start analyzing the role of vector fields in 
cosmology without the ordinary matter, and include the matter into consideration after the 
main features of vector fields are clarified. 

According to observations the Universe expands, and its large scale structure remains 
homogeneous and isotropic. Consider the space-time with the structure T 1 i? a!o and the 
metric 

do 

ds 2 = gjKdx'dx* = (dx ) 2 - e 2F ^ J^dx 7 ) 2 ( 7 ) 

depending on only one time-like coordinate x° = ct[l9]. Here d = 3 is the dimension of 
space. However, the derivations below are applicable for arbitrary do > 1. The metric tensor 
giK is diagonal. The uniform and isotropic expansion is characterized by the single metric 
function F(x°), and the rate of expansion is J™ = F' . The Ricci tensor is also diagonal: 

Roo = -d (F' 2 + F") J (8) 

Rjj = e 2F (F" + d F' 2 ), J>0. (9) 



The detailed analysis of vector fields in the cosmological metric ([7]) |16|] showed that the 
simplest case 

a^O, b = c = (10) 



is the most interesting one from the point of view of the Universe evolution. 

Depending on the sign of the invariant ^(pi the vector (ft 1 is either time-like {(f) 1 (f>i > 0), 
or space-like {(f) 1 (f>i < 0). In general relativity, while all coordinates are formally equivalent, 
one can choose the appropriate coordinate system where (f) 1 = either for I > 0, or for 
1 = 0. But it can not be done if the coordinate system is already chosen in accordance with 
some other reasons. In the cosmological metric (J7|) the coordinate x° is already specified, 
and whatever the sign of the scalar 7 0/ is, we have to consider 0/ having both space and 
time components. 

However, all space coordinates in the metric (j7|) are equivalent, and we can choose the 
coordinate x 1 along the space direction of the vector field. Then the vector 0/ has only two 
nonzero components 4>q, and (f>\. All other space components of the vector (f>j are zeros: 

(f)j = 0, I>1. (11) 

In the case (fTUj) the vector field equations @ reduce to 

W = -V'h, 1 = 0,1. 02) 

*^ a (^ + d0 ^ - e -^). 03) 
The energy-momentum tensor ((Sj) reduces to 

T IK = giK ($7« + V) + 2V (4>i4> K - 9iK(f)L(f) L ) ■ (14) 



As usual, the scalar 

V'{0)/a = m 2 

in the case (TTDT) can be designated as the square of mass of a vector field. 
A. Massless field 

The equations ( 1T2|) for a massless field, m = 0, are simply 



dx 1 

The divergence of the vector field (f> L L is a constant scalar 



0. (15) 



S L L = m = 0. (16) 
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The energy- momentum tensor (JUJ), 

T IK = g IK ($ 2 /a + V ), V' = 0, (17) 

acts in the Einstein equations (jBJ) as a simple addition to the cosmological constant: 

Rjk - \qikR + AgiK = 0, A = A - x ($jj/a + V ) . (18) 

Here Vo is the value of the constant potential V (0l0 l ) in the case of massless field (V = 0). 
The contribution of the zero-mass field to the curvature of space-time remains constant in 
the process of the Universe evolution. 
The metric 

ds 2 = (dxy~e ± ^^ ix °~ X " ) J2(dx 1 ) 2 , d >l (19) 

i=i 

is the self-consistent regular solution of the Einstein equations (JUJ), provided that 

A < 0. (20) 

F(x°) is a linear function; Xq is a constant of integration. The metric ffl9|) is called de Sitter 
(or anti de Sitter, depending on the sign definition of the Ricci tensor). It describes either 
expansion (sign +), or contraction (sign — ) of the Universe at a constant rate. In the case 
of sign + the rate of expansion 

H ^-lJ^T) (21) 

is called Hubble constant. In our 3-dimensional space 



H=d~A, d = 3. 
In general relativity the requirement of regularity f l20|) should replace the artificially 



imposed restriction (J3J) w 
negative energy problem 



rich people had been using for a long time in order to avoid the 



The zero mass vector field determines the constant rate of expansion. Available today 
properties of the so called dark energy (presently unknown form of matter providing the 
major contribution to the uniform isotropic expansion of the Universe) can be described 
macroscopically by the zero-mass vector field with a simple Lagrangian 

L = a (<P%) 2 - V . (22) 
7 



As long as the physical nature of vacuum is not known, the "geometrical" origin of the 
cosmological constant A and the "material" contribution to A by the zero-mass vector field 
can not be separated from one another. The combined action of the massless field and/or 
the cosmological constant is determined by the single parameter - Hubble constant ( )2ip . 



B. Massive field 



There is a principle difference between the massless and the massive longitudinal vector 
field. The massive field vanishes with time in the process of expansion, while the massless 
one does not [16]. If V 7^ the field equations (TT2"1) allow to express O an d 4>\ via $ : 

1 <9$ 



0,1. 



(23) 



V dx 1 ' 

Over the scales much larger than the distances between the galaxies the Universe is uniform 
and isotropic. Hence, there is no dependence on the space coordinates, and in accordance 
with (123]) <pi = 0. In the scale of the Universe the massive vector field is longitudinal: the only 
nonzero component 0o is directed along and depends upon the same time-like coordinate x°. 
The energy- momentum tensor ( Tl4l for the massive field reduces to [l6] 



T, 



IK 



gm (0o + d F'(f)o) 2 + (2S io 5ko ~ 9ik) m 2 ^ . 



dV 



Here the prime denotes the derivative d/dx° (except that V — 
massless vector fields can be of different physical nature. In the Einstein equations 



(24) 



Massive and 



-d (d - l)F' 2 + A = xa 
2 



% + d o F'0 o ) 2 



+ m 2 4>l 



{do - 1) F" + ^d (d - 1) F' 2 + A = xa 



b'o + d o F'0 o ) 



2 i2 

m O 



(25) 



(26) 



the massless field is taken into account via the constant A, and the massive one is described 



by the function (f>o 16| . Their applicability is restricted by the condition that the second 
and higher derivatives of the potential V (0l0 l ) can be ignored. Regular solutions of these 
equations are possible if A and a are of the same sign. It follows from (1251) . provided that 
F' can change the sign. 

The vector field equations 01]) reduce to the only one equation 



b'o + d F'4>o)' + m 2 4> = 0, 



(27) 



which is the consequence of the Einstein equations ( |25ll26|) due to the Bianchi identities. 
Extracting f )25|) from (126]) . we have 

F" = -Am 2 i (28) 
a — i 



The fourth order set of equations (12711281) with the boundary condition 

0o 2 + m 2 (pl = — at F' = 0, (29) 
xa 

following from (125]) . has the same solutions as the third order set ( 1251126]) . In the regular 
expanding solutions at x° — > oo , in accordance with (120]) . A < 0. Hence, a is also negative, 
and F" (128]) is positive: the massive longitudinal vector field makes the expansion of the 
Universe accelerated. Without ordinary matter the second derivative (128]) does not change 
sign. x° is a cyclic coordinate, and it is convenient to set F' = at the origin x° = 0. In 
view of the x° — > —x° invariance of the equations, it is clear, that the rate of expansion 
F' (x°) is a monotonically growing function between its limiting values F' {—oo) = —H in 
the past and F' (+oo) = H in future. The scale factor R = e F decreases with time while 
F' < 0, reaches its minimum, and grows when F' becomes positive. 

One of the two constants 0o and 4>' at x° = remains arbitrary within the boundary 
condition ( 129]) . In the case O (0) = the field 0o {%°) is a symmetric function, and in the 
case 0o (0) = it is an antisymmetric one. In these both cases F' (x°) is antisymmetric. If 
both constants </>o and <j)' at x° = are not zeroes, the regular solutions still exist, but there 
is no symmetry with respect to x° — > —x°. 

It is convenient to perform the further analytical and numerical analysis in dimensionless 
variables <t> and z: 

2UX T O , z = Hx°. (30) 



The equations (1271128 



d {d - 1) 

+ d o i r , 2 0) i2 + /i 2 = O, (31) 
F z , z = d fi 2 <p 2 (32) 



contain only one dimensionless parameter 



„ = |. (33) 
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(In dimensional units \i = ?§r)- Subscript |2 denotes d/dz. The boundary relations ( )29|) in 
terms of (0) , </> )Z (0) , and F a (0) are 

2 ,(O)+/i 2 2 (O) = l, F,(0) = 0. (34) 

F (z) enters the equations (j31|32p only via the derivatives, but not directly. For this reason 
Fq = F (0) remains arbitrary. The constant Fq only determines the scale of the space 
coordinates, and it does not change the structure of the metric and vector field. 



1. Asymptotic behavior at large x° 

At late times the back reaction of the massive field on the metric becomes negligible. In 
accordance with ( 132"]) F )Z)Z —> 0, and F jZ — > 1 at x° — > oo. The late temporal evolution of the 
massive field is described by (13T]) with F >z = 1. Its solution 



do , / f do^ 



<P(z) = C + e x + z + C^-\ \ ± = -JL±J\JLj -f, z^oo (35) 

is a linear combination of two functions, vanishing at z — >■ oo. The functions are monotonic, 
if /i < y, or oscillating with a decreasing magnitude, if /i > ^. If /i is small, the field 
decreases very slowly: 



(z) = C + exp ( — -^z ) , 2 — >■ oo, fj,<^ d . (36) 



2// 
"o 

For large and small \x the equations (13111321) can be solved analytically. 



2. Large n 

In the case of large //, 



H » 1, (37) 



the field (z) is a rapidly oscillating function as compared with F (z) . We search the solution 
of the field equation (13T1) in the form 

4> (z) = —ip (z) COS (flZ + <y?) . 
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In accordance with (I34j) ip (0) = 1, and the phase tp depends on the relation between the 
boundary values (0) and )Z (0) . The main terms ~ fi in fl3T|) disappear, the next order 
terms ~ 1 are: 

2</>, 2 + = 0. 

The remaining terms ~ <C 1 can be neglected. The solution is 

( Z ) = _ e ~MF-Fo)/2 cog / + s fi ^ 1 

Here F = F (0) . 

Averaging (132]) over the z-interval much larger than the period of vibrations fi' 1 , we get 
the following equation for F (z) : 

The solution is 

F (z) = F + 1 In cosh , F 2 (z) = tanh fe) , (38) 

, , . 1 cos (uz + ip) , . 

*^ = T^m- " >>1 - (39) 

Oscillations of (z) at large fi initiate weak vibrations of the function F (z) around 
the averaged value F (x) , see Figure HJ Red curve is the numerical solution for fi — 5, 
0,2 (0) = 0. Blue dashed line - analytical solution F z (z) (1381) . In Figures 2-left and 2- 
right (z) , found numerically for fi = 5 , practically coincide with found analytically (J39]) . 
The boundary condition >2 (0) = for the symmetric solution in Figure 2-left corresponds 
to if = in (131?1) . The antisymmetric solution in Figure 2-right (the boundary condition 
(J) (0) = 0) coincides with (pj) at <p = -tt/2. 

5. Case /x is very small 

If in the case /i<lwe neglect the term /^ 2 in the vector field equation ( 13TI) . then the 
vector field (z) is expressed via F (z) as follows: 

(*) = ^ (0) J* e do ^- Fo Uz 1 + (0)^ e - do[F ^~ Fo \ fi < 1. (40) 

The constants (0) and i2 (0) obey the boundary relation ( 1341) . 
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FIG. 1: Red curve is the numerical solution F )Z (z) for [i = 5, and F jZ (z) = tanh (^zj — blue 
dashed line; do = 3. 
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FIG. 2: Symmetric (left) and antisymmetric (right) <f> (z), found numerically for fj, = 5 and do = 3, 
coincide with found analytically (|39p . 



Symmetric and antisymmetric configurations 

In the symmetric case j2 (0) = (1301) reduces to 



-l e -d [F W -F ] j /i<<1 



(41) 
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FIG. 3: Symmetric solution F (z) — Fq, F tZ (z) , ficj) (z) found analytically (ji2]) for fj, -C 1. do = 3 

0(0) = /i" 1 in accordance with the boundary relation ( )34|) . Neglecting the term /z 2 (/> in 
the vector field equation ( 13T]) results in C_ = in the asymptotic behavior ( 13"5j) of (2) 
at z 00. On the contrary, omitting fi 2 <f) in (13"T1) in the antisymmetric case </> (0) = 
corresponds to C + = 0. The integral in ( HOI converges at z — )■ 00, and both expressions, ( HOI 
and ( 1331) . give the same result: (2) — > const at — > and z — >■ 00. Antisymmetric (2), 
acting as a massless field in the limit // — > 0, renormalizes the Hubble constant. Therefore 
it is reasonable to consider the contribution of the antisymmetric <fi (z) at fi — > as already 
included into A. 

In the symmetric case the metric is determined by the equation 

F )Z)g = d e~ 2d °™- F ^ 

and boundary conditions 

F )g (0) = 0, F (0) = F . 
The analytical solution for the symmetric configuration 

F 0) = F + -j- In cosh (do*) , F, (z) = tanh (d z) , (f) (z) = } , - . , /x < 1, (42) 

ct /i cosh [d z) 

describes the transition from the compression to the expansion, see Figure 3. 

With no ordinary matter the time interval of transition is of the order of Hubble time 
~ 1/ d H. At fi <C 1 it does not depend on the mass m of the massive field. The scale factor 
is 

R ( Z ) = e F{z) = e F ° [cosh (d z)} 1/d ° . 



13 




-WJ -29 



FIG. 4: Antisymmetric solution F (z) — Fq, F ;Z (z) , (j> (z) , found numerically for fi = 0.25, do = 3 
Without the ordinary matter the acceleration F )Z)Z is positive: 



cosh 2 (d z) 



> 0, 



Like an elastic spring, the symmetric longitudinal vector field enables the transition from 
compression to expansion. The kinetic energy of contraction completely converts at x° = 
into potential energy of the compressed vector field, and at x° > the energy is being 
released back in the form of the kinetic energy of expansion. 

An example of the antisymmetric solution (with the boundary condition (0) = and 
fj, = 0.25), found numerically for d = 3, is presented in Figure 4. 

Numerical analysis confirms, that at /i 1 the antisymmetric (z) decreases very slowly 
as at z — > oo, see ( )36|) and Figure 4 (right). The symmetric field (Figure 3 right) vanishes 
quickly, ~ exp (—doz). 



IV. EVOLUTION DETERMINED BY VECTOR FIELDS AND ORDINARY MAT- 
TER 

The above analysis of general properties of the equations (I31|32p facilitates clarifying the 
solutions of the Einstein equations with the ordinary matter taken into account. 



A. Dust matter approximation 



Applying the general relativity to the Universe as a whole it is natural to consider the 
ordinary matter (stars, galaxies, ...) as separated noninteracting subsystems located far from 

14 



one another. Averaged over the distances larger than the distance between the objects, the 
ordinary matter can be considered macroscopically as a uniformly distributed dust. As far 
as the ordinary matter does not violate the homogeneity and isotropy of the large scale 
structure of the space, all the components of its energy-momentum tensor, except T 00 = p, 
are zeros: 

Tik = P<5/o<W (43) 

In the process of expansion the averaged energy density of matter p depends only on x°. 
Tjk is a symmetric tensor. Due to the Bianchi identities its covariant divergence is zero, 



rnK _ 1 d / r-rpK\ ^ &9KL rpKL _ n 

h ' K ~ ~^gdx* W9lj ] ~ 2~d^ 1 ~ U - 
In the cosmological metric (j7j) the covariant divergence QHj) of the tensor (|43p is 



(44) 



T K _ -d F V r) r 

Thus, pe d ° F = const. If we accept that p (x *) = po is the averaged density of the ordinary 
matter now, then 

p (x°) = poe-*^* ) , 
and the present moment x°* is defined by 



F (x *) = 0. 



(45) 



B. Einstein equations 



With the dust matter taken into account, the Einstein equations ( l25|26|) are: 

2 



1 



d (d — 1) F + A = na 



dx° 



°+d F'<p ) +mVo 



+ xp e 



(d - 1) F" + -d (d - 1) F' 2 + A = xa 



-d F 



2 j. 2 

m <p, 



In dimensionless variables (130 



Fl-1 



±F + f 2 - 1 

do ' z > z > z 



+ d F z <p) 2 - p 2 <p 2 + Qe~ doF , 
-{<p, z + doF^f + p 2 4> 2 . 
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(46) 
(47) 



The parameter Q, 

_ _xpo _ 2xp 

A ~d (d -l)H^ 

denotes the ratio of the energy density of the ordinary matter to the density of the kinetic 
energy of expansion. The vector field equation (I3TI) remains the same, and instead of the 
equation ( |32i) we now have 

F,, = - ^e~ d ° F . (49) 

This equation resembles the Newton's law: acceleration F )ZjZ is proportional to the "repulsing 
force" do/j, 2 (j) 2 minus the "attracting force" ^ L e~ d ° F . Altogether the regular solutions of the 
set f )3T|49|) with the boundary conditions 

(0) + /i 2 2 (0) = l + Qe- doFo , F fg (0) = 0, F(0) = F , A < 0, (50) 

contain five dimensionless parameters: /i, fl, F , <p (0) , and <p >z (0) . In view of the bounding 
relation ( 150]) . four of them remain independent. 



C. Regular cosmological solutions 

The equations ( l3Tf49l with the boundary conditions ( 1501) are easily integrated numerically. 
The regular solutions are free from fine tuning. Moreover, the existing parametric freedom 
results in a great variety of possible configurations. Among the four independent parameters, 
fl (jlSJ) can be estimated better than others. Substituting do = 3, x = 6.67 x 10~ 8 cm 3 /g 



sec 2 , H = 2 x 10~ 18 sec -1 , po — 2 x 10 -31 g/cm 3 (universal density of luminous matter 20]) 
we have Q ~ 10 -3 . Numerical analysis shows, that if both <fr (0) , and <fi tZ (0) are not zeros 
at a turning point F yZ = 0, then there can be other even more sharp turning points with <fr jZ 
more close to zero. With fixed values of two parameters, (fl = 10~ 3 , jZ (0) = 0) , the role 
of \x at fixed F = —3 in the process of evolution is demonstrated in Figure 5, and the role 
of F at fixed \x = 3 is shown in Figure 6. 

The rate of evolution F yZ (z) for fi — 0,1, and 2 is presented in Figure 5 by brown, blue, 
and red curves, respectively. In all symmetric configurations the global transition from 
contraction to expansion takes place at z = 0. At fixed Fq = —3 after the transition the 
rate of expansion reaches its first maximum F 2 (z max ) > 1 at some z max > 0. At z = z m3iX 
the acceleration turns into deceleration. If /i = (brown curve in Figure 5) F z (z > z m£LX ) 
decreases monotonically to F )g — 1 as z — >■ oo. If /i > a next point of extremum z min > z max 
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FIG. 5: The rate of expansion F )Z (z) . [i = 0, 1, 2 - brown, blue, red curves, do = 3, O. = 0.001, 
F = -3. 



appears, where deceleration turns back into acceleration. The bigger is /i, the deeper is the 
minimum of F jZ at z = z min , - compare the blue (/i = 1) and red (fj, = 2) curves in Figure 
5. The number of subsequent maxima and minima grows with increasing \i. 

Parameter F determines the maximal speed of the outburst F )Z (z max ) in the contraction- 
to-expansion transition. The rate of expansion F 2 (z) for F = —3, —4, and —5 at fixed \i — 3 
is presented in Figure 6. The maximum F iZ (-2 m ax) grows exponentially with the increasing 
negative value of Fq. It resembles inflation, except that there is no singularity. Some authors 
call the regular contraction-to-expansion transition "nonsingular bounce" 21|, 22], 23]. In 
the literature there are attempts to find a self-consistent model in order to explain from a 
unified viewpoint the inflation in the early Universe and the late-time accelerated expansion 



241 ]. However, one should keep in mind that the dust matter approximation is applicable 



until the galaxies are noninteracting systems, located at far distances from one another. 

In the case of small fi 1 the transition from contraction to expansion can be described 
analytically. Substituting ( 14T]) into ( 149|) . we exclude (f) and come to the single equation for 
F : 

d Q 



F z , z = doe- 8 *^ 

Its solution with the boundary conditions (1501) is 

1 



-d F 



F Z) 



F + — In 

do 



1 + ^ne~ doFo ^j cosh (faz) - hle 



d F 



n < 1. 



(51) 



(52) 
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FIG. 6: The rate of expansion F tZ (z). fx = 3. Fq = —3, —4, —5 - brown, blue, red curves, do = 3, 

n = o.ooi 
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FIG. 7: Variation of F z (z) in the vicinity of the turning point, found numerically for /i = 10 and 
Fq = —10, coincides with ([53]) 



For the rate of evolution F z (z) and for the scale factor R (z) 



we get 



F z (z) 



sinh (doz) 



cosh (d z) - (1 + n edoi?0 ) 



-i • 



3 doFo + In ) cosh (d 2) - in 



(53) 



(54) 



2 J y ' 2 

Analytical solutions f )52ll54|) . derived for /i C 1, are as well applicable in the vicinity of the 
transition for fj, ~ 1 if |Fo| ^> 1. See Figure 7, where the variation of F jZ (z) in the vicinity 
of the turning point, found numerically for fj, = 10 and Fq = —10, coincides with (l5"3l . It is 
because for very large negative F the width of the contraction-to-expansion transition Az 
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FIG. 8: The vector field (z), the rate of evolution F tZ (z) , and the metric function F (z), found 
numerically for d = 3, <f> fZ (0) =0,n= 10, F = -4, and O = 0.001 

is very narrow: 

Az 2 — e doFo/2 , F < 0, |F | > 1. 

According to the recent analysis of the Hubble space telescope data, the expansion of 
the Universe switched from deceleration to acceleration at about the half of the age of 
the Universe [25j. In the analytical solution ( 152]) F iZ>z (z) is negative at z > z max : if /i = 
0, the expansion goes with deceleration. The switching from deceleration to acceleration 
(appearance of the minimum of F jZ (z) at z — z m i Q , see Figures 5,6) is a hint that the 
symmetric longitudinal vector field is massive, /i > 0. At [i > y the field (z) becomes an 
oscillating function, and the number of minimums of F >z (z) grows with growing /i. The vector 
field <f){z), the rate of evolution F >z (z) , and the metric function F(z), found numerically 
for (f> tZ (0) = 0, do — 3, [A — 10, F = —4, and Q = 0.001, are shown in the Figure 8. For this 
set of parameters, in accordance with ( 145| . the today's "date" is z* = Hx* = 3.45. 

The damping oscillations of the vector field <fi (z) give rise to the oscillations of F z (z) and 
result in slight variations of F (z) in the vicinity of the compression-to-expansion transition. 
Practically, for /i 3> 1 the oscillations vanish at z > 1. 

Another possible origin of the deceleration-to-acceleration switching at about the half of 
the age of the Universe (within our approach) can be a slight antisymmetric contribution of 
0,2 (0) 7^ in the boundary conditions ( |50l) . Numerical analysis shows, that a small nonzero 
0,z (0) , |0,2 (0)| <C |0 (0)| , can also lead to the appearance of a minimum. 
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D. Regular solutions with positive A 



There is an important difference between the boundary conditions (j34p and f )50|) . The 
relation (1341) can be satisfied only if A < 0, provided that a < 0. Appearance of the term 
Q e -d F j n ^ ac [ m i^ s j^e solutions with positive A. If A changes sign, then H f[2"Tj) be- 
comes imaginary. The equations ( 13T1H91) are invariant against H — > iH, but the boundary 
conditions fl50l) are not: 



(55) 



2 z (0) + /i 2 2 (0) = -l + Qe- doFo , F z (0) = 0, F(0) = F o , A > 0. 



The necessary condition for regular solutions with A > is the existence of an extremum 
moment (F x (0) = 0) with the energy density of the ordinary matter exceeding the kinetic 
energy of expansion: 

-d F _ (0) 



A 



> 1, F z (0) = 0, A > 0. 



In the case A > 0, /i <C 1 the analytical solution of the equations (I31|49p with the boundary 
conditions fl55|) is 



F (z) = F + 1 In 
d 



1 - ^ne~ doFo J cos (doz) + -Oe^ 



<f)(z) = fi 1 

The scale factor fz) and the rate of evolution F z (z 



1 - ~n e - doFo j cos (d„z) + ^fie- doFo 



-fie - *^ - ( ifie- doFo - 1 ) cos (do^) 



F z (z) 



sin (do-^) 



(l_2 e d Fo) ^cos^qz)' 

are periodical functions with no singularities, see red curves in Figure 9. For the values of 
the parameters do = 3, p. — 0.02, Fo = —3, Qe~ d ° F ° = 1.032 there is no difference between 
the curves found numerically and analytically. 

Without the massive field (0 = 0) the solutions with positive A are possible only if the 
parameters are fine tuned ( xp (0) = A ) : 



F(z) 



F + — In cos — , 
do 2 
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R(z)/R(0) 
1 




FIG. 9: Left: scale factor R(z)/R(0). Right: rate of evolution F tZ (z). Red curves - numerical 
(coinciding with analytical) solutions for do = 3, fj, = 0.02, Fq = —3, Qe~ d ° F ° = 1.032; blue curves 

- (E7D and (Ea, n e - d ° Fo = 1. 



F z (z) 



R(z) 



tan 



d z 



f #- d o z 
e u cos d o . 



(56) 
(57) 



The scale factor R (z) ( |57|) and the rate of expansion F z (z) (|56|) are presented in Figure 



9 (blue curves). These "fine tuned" ^e Fo = ^ 



do(d -l)H 2 



1/do 



singular solutions have peri- 



odical singularities at z = z n = ^ (1 + n) , n — ±1, ±2, ... . In the vicinity of each singular 
point z n = Hx° n1 as well as at H — > 0, the Hubble constant H drops out, and the scale factor 
( 1571) (in the ordinary units x° = z/H) reduces to 



R(x°) 



d Hp 



l/do 



X 



„0\ 2 / d o 



X 



< 1. 



(58) 



n 



2(d Q -l) 

At <i = 3 (158|) reproduces the scale factor of the Friedman-Robertson- Walker [2j cosmology 
with dust matter in the plane space geometry. The longitudinal vector field <fi ^ removes 
the singularities, see red curves in the Figure 9. 



26| as an alternative 



The idea of "oscillating Universe" is actively supported by Lessner 
approach to cosmology. Meanwhile, I would rather not call it alternative, for all the above 
solutions, including the oscillating ones, are derived completely within the frames of the 
standard Lagrange approach and General relativity. 
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V. SUMMARY 



Big Bang singularity is not the inevitable property of the Universe evolution. It is a 
consequence of our not knowing the physical nature of the dark sector, including the origin 
of its ability to resist compression. Nevertheless, there is a possibility for macroscopic 
description of the Universe evolution within the frames of the Einstein's theory of general 
relativity. The simplest non-gauge vector field with the Lagrangian 

L = a {[ct> M M ) 2 - m 2 <^) - V (59) 

can be that missing link in the chain, necessary to understand the mechanism of accelerated 
expansion of the Universe, and avoid, better say - resolve, the Big Bang singularity. 

It is rather involuntarily, but the modern interpretations of the observational data are 
based on the idea of the Big Bang birth of the Universe. The cosmic background radiation, 
among other phenomena, definitely testifies that the Universe had been strongly compressed 
in the past. But how strongly? That is the question. The information from the past, coming 
to us with electromagnetic waves, tells us only about the phenomena that happened after 
the Universe became transparent. The far extrapolation to the Plank's era is based on 
the assumption that the singularity is the inevitable property of cosmological solutions of 
the Einstein equations. It is so for the solutions, taking into account only the ordinary 
matter (including electromagnetic radiation). The discovery of the accelerated expansion 
strictly pointed on the existence of hidden sector, able to resist the compression. Though 
the presented above solutions are regular, they do not restrict the value of compression 
from above. The degree of maximum compression is determined by a free parameter F . At 
large negative F the regular transition from compression to expansion looks like inflation. 
However, one should keep in mind that the dust matter approximation is applicable if the 
galaxies are located at far distances from one another. 

The macroscopic theory can help understanding some specific features of the Universe 
evolution, but it does not determine the physical nature of what we call "dark sector". 
The vector fields can be considered as appropriate tools for macroscopic description of the 
Universe evolution. From my point of view, it is more important what the fields do, than 
how we call them. The previous analysis confirms that the zero-mass vector field is adequate 



for the description of dark energy 16J . In some papers the dark energy is considered as the 



origin of acceleration. However, zero-mass fields can be responsible for contractions and 
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expansions only at a constant rate. The acceleration is connected with massive vector fields. 
Being massive, they can not be associated with dark energy. At the same time, their energy- 
momentum tensor differs from the one of the ordinary matter. Therefore, the massive fields 
should correspond to dark matter (at least, partly). 

From the General relativity viewpoint, all three kinetic terms in the Lagrangian ([2]) have 
equal rights. In flat spacetime the antisymmetric term ~ (b — c) corresponds to electro- 
magnetic field (photons). The symmetric term ~ (b + c) has attention to the gauge vector 
particles There is no reason why the term ~ a should be "more equal than others". So, 
the Lagrangian f l5^|) also deserves to be associated with some particles, existing in nature. 
In accordance with the Subsection "Regular cosmological solutions", the observed [25j point 
of minimum z min , where the deceleration turns back to acceleration (see Figures 5,6), could 
correspond to \i ~ 1. It is worth trying to detect an extremely light particle with the rest 
energy mc 2 ~ hH ~ 1(T 33 eV. The trouble is that the massive vector field displays itself via 
gravitation (it curves the spacetime), but there is no evidence of its direct interaction with 
the ordinary matter. 
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